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Abstract 

We construct a discrete analogue of the integrable two-dimensional Dirac operator and describe 
the spectral properties of its eigenfunctions. We construct an integrable discrete analogue of the 
modified Novikov-Veselov hierarchy. We derive the first two equations of the hierarchy and give 
explicit formulas for the eigenfunctions in terms of the theta-functions of the associated spectral 
curve. 



^ 1 Introduction 

The purpose of this paper is to construct a discrete analogue of the modified Novikov-Veselov hierarchy 
I— land its algebro-geometric solutions, and to describe the spectral theory of the corresponding discrete 
I Dirac operator. 

>■ ■ The modified Novikov-Veselov (mNV) hierarchy is an integrable hierarchy of equations introduced 
^ by Bogdanov in [T] , [2] as a special reduction of the Davey-Stewardson equation. The equations of 
the hierarchy have the form of Manakov L, A, i?-triples 

cn ; 

^ = [L,A^]-B^L, (1.1) 
^ ifvhere L = D is the two-dimensional Dirac operator 

I ^''-(^ ')(::)• <-) 

and An and Bn are (2 x 2)-matrix differential operators. The mNV hierarchy describes deformations 
of the Dirac operator that preserve the zero energy level, i.e. isospectral deformations of the equation 

^^ = 0. (1.3) 

The first equation of the hierarchy has the form 

Ut = (^U;,^^ + 3U^V + ^UV^^ + (^U^2z + 3UzV + ^UV2^ , = (u^)^. (1.4) 

In [3], [H Taimanov constructed algebro-geometric solutions of the mNV hierarchy and described 
the spectral theory of the Dirac operator (11. 2p . In recent times, the mNV hierarchy and its algebro- 
geometric solutions have attracted significant attention due to their applications to the classical theory 
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of two-dimensional surfaces in three-dimensional Euclidean space, and in particular to the Willmore 
conjecture (see the survey [5j for an extensive bibliography). 

It is possible to consider a more general two-dimensional Dirac operator of the form 



D 



u d 
—d V 



) 



(1.5) 



The spectral theory of the two-dimensional Dirac operator (11.51) is equivalent to that of the two- 
dimensional scalar Schrodinger operator in a magnetic field 



The reduction of the Dirac operator (I1.5P to the form (II. 2p corresponds to a reduction on the 
Schrodinger operator in which the functions U and V satisty the relation 



The analytic properties of Baker-Akhiezer functions which describe general Schrodinger operators of 
the form (II. 6p that are integrable on the zero energy level were formulated in [lOj. The reductions 
on the algebro-geometric data that describe the potential Schrodinger operator {V = 0), which is the 
auxiliary operator for the Novikov-Veselov hierarchy, were found in [8], [9]. 

The problem of constructing an integrable discretization of an integrable differential equation is 
not well-posed and does not have a universal solution. However, there are several methods in soliton 
theory that allow us to construct integrable discretizations. Most of them are based on constucting a 
discrete analogue of the auxiliary linear problems, which involves an appropriate deformation of the 
analytic properties of the solutions of these linear problems. 

In the finite-gap case, the eigenfunction of the auxiliary linear differential operator, known as 
the Baker-Akhiezer function, is defined on an algebraic Riemann surface and is required to have 
exponential singularities controlled by the continuous variables at one or more marked points of the 
surface. To construct a discrete analogue of the operator, we replace each exponential singularity 
with a pair of meromorphic singularities consising of a pole and a zero of the same order, which 
we consider as the discrete variable. This deformed eigenfunction then satisfies a infinite system of 
linear difference and differential equations, whose compatibility conditions are the discretization of the 
original integrable hierarchy. This method was used for constructing algebro-geometric solutions of 
the Ablowitz-Ladik equation [11], [12], which is a discretization of the nonlinear Schrodinger equation, 
and for constructing Darboux-Egoroff lattices, which are the discrete analogue of Darboux-Egoroff 
metrics [S]- 

Using this approach, Grushevsky and Krichever have given an algebro-geometric construction of an 
integrable discretization of the two-dimensional Schrodinger operator (II. 6p . In the second paragraph, 
we describe a matrix variant of this construction, which leads to a two-dimensional matrix difference 
operator of the form 



where the functions ipi and the coefficients of the operator are functions of two discrete variables 
n,m & and Ti, T2 denote the translation operators in the discrete variables. The operator D, which 
we call the discrete Dirac operator, can be considered as a discrete analogue of the general Dirac 
operator of the form (II. 5p . 



H = dd + Vd + U. 




V = -dlnU. 



(1.7) 
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The coefficients of a discrete Dirac operator (11. 8p depend, up to gauge transformation, on two 
arbitrary functions of the discrete variables. In the second paragraph, we show that a discretization 
of the algebro-geometric data corresponding to operators of the form (II. 2p leads to operators whose 
coefficients depend on only one arbitrary function, namely operators of the form 



Dtp 



T2 
Ti 



a P 
(3 a 



^1 



where the coefficients satisfy the relation 



a 



P' 



;i.9) 



;i.io) 



In the third paragraph we introduce time dependence into the eigenf unctions and construct an inte- 
grable hierarchy of isospectral deformations of the zero energy level of the operator (11.91) . We call 
this hierarchy, which has the form of Manakov L, A, 5-triples, the discrete modified Novikov-Veselov 
hierarchy. 

In the fourth paragraph we derive the explicit form of the first two equations of the hierarchy 
(equations (14.221) . (I4.24p . (14.261) ). The ffist equation has the following form: 



^/^g2(/3(n— l,m+l) g2(/p{n— l,m) ^ ^Q—2ip(n,rri) g— 2i^(n,m+l) ^ 



where the two functions satisfy the non-local relation 

(p{n, m + 1) — ^{n, m) = tp{n + 1, m) — ip{n, m). 



1.11) 



1.12) 



In the final paragraph we give explicit formulas for the Baker-Akhiezer functions in terms of theta- 
functions associated to the spectral curve. 



2 Reduction of general discrete Dirac operators 

Consider the following discrete linear equation 



T2 
Ti 



a P 
7 6 



^1 

i'2 



0, 



(2.1) 



where ip = {ipiin, m) , ilj2{'n, m))'^ is a vector function of two discrete variables n,m ^Z, and 



a P 
7 6 



a{n,m) P{n,m) 
7(n, m) 6{n,m) 



(2.2) 



is a (2 X 2)-matrix function of the discrete variables. We call D the discrete Dirac operator. We use 
Ti and T2 to denote the translation operators in the discrete variables 



Ti/(ra, m) = f{n + 1, m), T2/(ra, m) = f{n, m + 1), 



(2.3) 



while ti and ^2 will be used to denote the translated functions, so that for example Ti{fg) = {tif){tig). 
In this chapter, we construct algebro-geometric solutions of equation (12. ip and some of its reductions. 

The main method of constructing algebro-geometric solutions of linear differential or difference 
equations such as (12. ip is to consider functions ifji defined on an auxiliary Riemann surface, called 
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the spectral curve, and having certain prescribed singularities on that curve. Generally, to construct 
solutions of difference equations, we consider functions that are meromorphic on the spectral curve 
with prescribed pole singularities, while constructing solutions of differential equations requires us to 
consider functions with prescribed essential singularities, called Baker- Akhiezer functions. 

Let X be a smooth Riemann surface of genus g. We consider the following data on X: 
Data A. 

• Four distinct marked points P^, on X. 

• Local parameters zf = (kf)~^ defined in some neighborhoods of these points. 

• An effective divisor D = ■ji + ■ ■ ■ + 7g_|_i of degree g + 1 on X, supported away from the marked 
points, which satisfies the following condition of general position: 

h\D + {n- 1)P+ - nP{^ + (m - 1)P^ - mP^) = for all n,meZ. (2.4) 

To construct solutions of equation (12.11) . we consider spaces of meromorphic functions on X with 
singularities controlled by the discrete variables: 

^n,ni = H°{D + nP+ - nP{ + mP^ - mP^) C Mer(A), n, m G Z. 

The Riemann-Roch theorem implies the following 

Proposition 1 Suppose that X is an algebraic curve with data A defined above. Then each of the 
spaces \l/n,m is two-dimensional: 

dim "^^,171 = h°{D + nP^ - nP{ + mP:^ - mP^) = 2 for all n,m eZ, 

the intersection of two of these spaces at adjacent lattice points is one- dimensional: 

dim n \&n,m-i = h°{D + nP^ - nP{ + (m - l)P^ - mP^) = 1 for all n,m eZ, 

dim n \l/n-i,m = h°{D + {n - 1)P^ - nP{ + mP:^ - mP^) = 1 for all n,m eZ, 
and these two one- dimensional subspaces o/\l/„.m span the entire space, i.e. their intersection is trivial: 

dim -^n,m n '^n,m-i H ^„-i,m = h° (D + {u - 1)P^ - uP^ + (m - l)P^ - mP^) = for all n,m eZ. 

Therefore, we can fix a basis ipiin, m, P),%l)2{n, m, P) in each of the spaces \l/ri,m by letting ipi{n, m, P) 
be any non-zero element of '^n,m H '^n,m-i, and letting ilj2{n,m, P) to be any non-zero element of 

ipi{n, m, P) e H^{D + nP^ - nP{^ + (m - 1)P+ - mP^) - {0}, (2.5) 

ij2{n, m, P) e H°{D + {n - 1)P+ - nPf + mP/ - mP^ ) - {0}. (2.6) 
The principal observation concerning these functions can be summarized in the following statement: 
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Proposition 2 Suppose that X is a Riemann surface with data A as defined above. Then there 
exist functions a{n,m) , (3{n,m), 'y{n,m), 6{n,m) such that the functions ipi{n,m, P) and ilj2{n,m, P) 
defined by ( l^.5]) -( fl7^) satisfy the Dirac equation: 



Dip 





^ 








-( 



a(n, m) f3{n,m) 



tpiin, m, P) 
ilj2{n,m,P) 



0. 



(2.7) 



Proof. Indeed, by construction, both ipi{n, m + 1, P) and ip2{n+l, m, P) actually lie in the space '^n,m, 
hence they can be expressed as linear combinations of the basis functions ipi{n, m, P) and ip2{n, m, P), 
which is equivalent to saying that the satisfy the Dirac equation ( ]2.7p . 

Therefore, a Riemann surface X together with the additional data given above allows us to construct 
a family of solutions {ipi{n, m, P),ilj2{'n, m, P))^ of the Dirac equation f l2.ip . parametrized by the points 
P of X. 

In order to construct reductions on the Dirac equation (12. 7p . we first express the coefficients a{n, m), 
P{n,m), '~f{n,m) and 6{n,m) in terms of the principal parts of the basis functions at the marked 
points. In terms of the chosen local coordinates, the basis functions ipi{n, m, P) and ip2{n, m, P) have 
the following expansions at the marked points, where k denotes the appropriate local parameter /c!|_: 



ipii^n, m, P) 



a^{n,m)k~''^ + 0{k-''-^ 

Oik""-^), 
02 {n, m)k~™- + 0{k 



m— 1^ 



ip2{n,m,P) 



0(A;"-i), 
6r(n,m)A;-" + 0(A;-"-i) 
bt{n,m)k"^ + olk 



as 


P - 




as 


P - 


-Pi 


as 


P - 




as 


P - 


-P2 


as 


P - 


-Pi 


as 


P - 




as 


P - 


.Pi 


as 


P - 


-.P2 



(2. 



(2.9) 



b2in,Tn)k~"' + 0{k- 

where the af{n,m) and bf(n,m) are functions of the discrete variables n and m. Considering the 
Dirac equation (12. 7p near the marked points P^,P^ gives us the following system of equations (in 
what follows, we usually suppress the indices n and m and replace them with the shift operators ti 
and ^2): 



t2aj 



aa 



1 5 



^1 ' 







(2.10) 





aa]" + Pb^ , tib2 
aa2 + ^0,2, tib2 

The functions ipi and ip2 have so far been defined up to multiplication by a constant factor dependent 
on n and m. We impose the following additional condition on the functions ipi and ip2'- 



7«i + 
6bt, 

1^2 +^^2 



btb2 



1. 



(2.11) 



It is easy to show using (I2.10p that these conditions imply the following relations on the coefficients 
a, f3, 7, 6: 

a 6 (t2af)(tib2^ 
5 " ~ " 



a6 — P'j 



±1. 



(2.12) 



a cti &2 

Condition (12. lip defines the constants af and 6^, and hence the functions tpi and 1^2, only up to a 
factor of ±1 that depends on n and m. This allows us to impose the following additional condition on 
the functions ipi and 4'2- 

{t2at){tib2) = atb2. (2.13) 
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In other words, we can choose the sign for the function ijj2 arbitrarily, and then choose the sign for 
the function ipi using the above relation. With this condition, the sign in equation fl2.12p is positive. 
Therefore, reductions ( 12. lip and (I2.13P impose the following relations on the coefficients of the Dirac 
operator (12. 7p : 

a5-p-f = l, a = 5 (2.14) 

In other words, the coefficients of a general Dirac operator of the form (12. 7p depend, up to gauge 
equivalence, on two arbitrary functions of the discrete variables. 

We now introduce a reduction under which the coefficients of the Dirac operator ( 12. 7p depend on 
only one function of the variables n, m. Suppose that, in addition to data A described above, the 
spectral curve X has the following 
Data B. 

• A holomorphic involution a:X X that interchanges the marked points and the local param- 
eters at the marked points as follows: 

a{Pt) = P^, o{kf) = kf. (2.15) 

• A meromorphic 1-form uj on X which has simple poles at the marked points with residues 
±1 and no other singularities, whose zero divisor is i5 + <y{D), and which is odd with respect to 
the involution. 

Consider the meromorphic 1-form ipiln, m, P)ilj2{n, m, a{P))uj{P). Comparing the singularities of the 
three terms, we see that this 1-form has simple poles at P^ and P2 with residues a^bi and — a^6^, 
respectively, and no other singularities. Hence, the existence of the additional data above implies that 
the coefficients of the functions tpi and 1^2 satisfy the following additional condition: 

a+fe- = a^6^ (2.16) 

Using (I2.10p and (12. lip , it is easy to show that this condition implies the following additional relation 
on the coefficients of the Dirac operator: 

/5 = 7. (2.17) 

Using the involution a we can rewrite the normalization conditions (12. lip and (I2.13P in the following 
equivalent form: 

V'i(m(a(p))ip^p+ = i, (2.18) 

^2(P)^2(^(P))|p=P+ = l, (2.19) 



t2MP) 



MP) 

Therefore, we can summarize the result of this reduction as follows. 



P=Pt UP) 



(2.20) 



Proposition 3 Suppose that X is a Riemann surface with data A and data B as defined above, and 
suppose the functions ipi{P) and ip2{P) defined by Ii2.5\) and i\2.6\) satisfy the normalization conditions 
i\2.18\) - f2.2U\) . Then there exist functions of the discrete variables a and (3 that satisfy the relation 

a^-(3^ = l (2.21) 

and such that the functions ipi{P) and ip2{P) satisfy the discrete Dirac equation: 



Dip 



T2 \ _ / a /5 
Ti i I /? a 



MP) 

MP) 



0. (2.22) 
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We now construct a further reduction of the discrete Dirac equation (12.221) which is the discrete 
analogue of the real- valued reduction in the differential case. Suppose that, in addition to data A and 
data B above, the spectral curve X has the following 
Data C. 

• An anti-holomorphic involution t:X ^ X that interchanges the marked points and acts on the 
local parameters at the marked points as follows: 

r(P±) = P±, r(P±) = P±, T{kt) = kl r{kt) = kf. (2.23) 

• A meromorphic function /(P) on X with divisor (/) = D — t{D) satisfying the conditions 

f{P)f{T{P)) = -1 for all P e X, f{P^)f{P{) = 1. (2.24) 

For a function f {71,171) of the discrete variables, we introduce the notation f*{n,m) = f (171,71). 
Consider the two functions V'2('^5 '^(■P)) ipi{7i, m, P)f{P). Both these functions are meromorphic 
and lie in the one-dimensional space H^{t{D) + {71 — 1)P2 ~ nP2 + ^Pi — rnPi), hence there exists 
a function C(n, m) of n and m such that 

Mm, 71, t{P)) = ipi{7i, 771, P)f{P)C{7i, m). (2.25) 

Considering this equation at P = P^ and P = P^ and using conditions (12. lip and ( 12.24p . we see that 

C(n, mf = 1 for all ra, m G Z. (2.26) 

We recall that the function 1^2 was normalized by condition (12. lip , which specifies it up to multiplica- 
tion by a factor ±1 dependent on n and m. Therefore, we can choose this factor in such a way that 
C{n,7n) = 1 for all n and m, in other words we may impose the additional following condition: 

V'2(m, n, r(P)) = ^n, m, P)/(P). (2.27) 

Equation ( 12.24p then implies that the functions ipi and ip2 chosen in this way satisfy the following 
relations: 

V'2(m,n,r(P)) =V^i(n,m,P)/(P), V^i(m, n, r(P)) = -V^2(n, m, P)/(P). (2.28) 

Plugging these relations into the reduced Dirac equation f l2.22p gives us the following relations on the 
coefficients of the operator: 

a* = a, 13* = -13. (2.29) 
We summarize the results of this reduction in the following proposition: 

Proposition 4 Suppose that X is an algebraic curve with data A, B and C as defined above, and 
suppose the functions ipi{P) and ip2{P) defined by Ili2.5\) and Ii2.6\) satisfy the normalization conditions 
I^2.18\} - f2.20\} and { 2.27\ ). Then there exist functions of the discrete variables a and (3 that satisfy the 



relations 

a^-P^ = 1, a* = a, (3* = -/3 (2.30) 
that the functions ipi{P) and ip2{P) satisfy the discrete Dirac equation: 



T2 \ _f a (3 
Ti i \ /3 a 



MP) 
MP) 



(2.31) 



7 



3 The discrete modified Novikov-Veselov hierarchy 



In the previous section, we constructed algebro-geometric solutions of the discrete Dirac operator 
(12. 7p and its reductions (12.221) and fl2.3ip by considering spaces of meromorphic functions \l/„,m, on an 
algebraic curve X with poles and zeroes determined by the numbers n and m. In this section, we 
embed these meromorphic solutions into a family of transcendental functions, called Baker- Akhiezer 
functions, and construct a hierarchy of commuting flows on the space of these functions. The set 
of compatibility conditions of these flows is the discrete analogue of the modified Novikov-Veselov 
hierarchy. 

Let r = {t}, rf, s = 1, 2, . . .} G C°° © C°° denote two sequences of complex numbers, only finitely 
many of which are non-zero, which we think of as continuous time variables. We construct deformations 
^n,m,T of the function spaces '^n,m constructed in Section 2 by considering functions which in addition 
have essential singularities at the marked points controlled by the times r. 

Proposition 5 Suppose that X is an algebraic curve with data A and data B given as in the previous 
section. Denote by X = X — — P{ — P2 — P2 the curve X with the marked points removed. 
Consider the space '^n,m,T ^ Mer{X) of functions on X defined by the following conditions 

1. For all ilj{n,m,T; P) G ^n,m,r we have + D >0, where (/) denotes the divisor of f . 

2. At the marked points P^ the elements ^p{n,m,T; P) of'^n,m,T have essential singularities of the 
following form, where by k we denote the appropriate local coordinate kf : 

^(n, m, r; P) = exp ^± ^ j 0(A;±") as P ^ Pf 
'4j{n,m,T]P) =exp (±^r^kA Oik^"") as P Pr 



± 



(3.1) 



s=l 



Then each of the spaces '^n,m,T is two-dimensional: 

dim '^n,m,T = 2 for all n,m E 7^, (3-2) 

the intersection of two of these spaces at adjacent lattice points is one- dimensional: 

dim '^n,m,T n '^n,m-i,T = 1 for all n,m E Z, (3.3) 

dim\E'„^m,r n '^n-i,m,T = 1 for all m G Z, (3.4) 
and these two one- dimensional subspaces of'^n,m,T span the entire space, i.e. intersection is trivial: 

dim '^n,m,T n \E'„,m-i,r H \E'„-i,m,r = for all n,m eZ. (3.5) 

Proof. The proof of this proposition is a standard application of the Riemann-Roch theorem. 

This proposition allows us to define functions ipi{n,m,T; P) and ip2{n,m,T; P) using the same 
relations as in Section 2. We observe the normalization conditions (I2.18I) - (I2.20I) can be applied to 
elements of '^n,m,T, since the exponential singularities cancel out. 
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Proposition 6 There exist unique functions ipi{n,m,T; P) and ip2{n,m,T; P) that form a basis for 
the vector space \l/„,m,T such that 

^l(n, m, r; P) G ^n,m,r H ^n,m-l,r " {0}, (3.6) 

z/'2(n, m, r; P) G ^n,m,r n ^„-i,™,, - {0}. (3.7) 

and which satisfy the normalization conditions l{2.18\) - l[2.20\) . These functions satisfy the discrete Dirac 
equation 



Dip 



T2 \ _f a (3 
Ti ) \ 13 a 



^1 
^2 



0, (3.^ 



where a and f3 are functions of the variables n, m, and r that satisfy the condition 

a^-f3^ = l. (3.9) 

In Section 6, we give explicit formulas for the functions ipi in terms of theta-functions. 

We now show that these functions satisfy a system of commuting linear equations. Let 91 denote 
the ring of functions in the variables n, m and r. We consider the ring D = 9^[Ti, T^"^, T2, T2~^] of finite 
difference operators with coefficients in 9^, and the ring 3Jl of (2 x 2) matrix operators with coefficients 
in O. By ip we denote the column vector [ipiln, m, r; P, 'ip2{'n, m, r; P))'^ ■ 

Proposition 7 There exist unique matrix difference operators Al in DJl 

a: = (<' . = 1.2. (3,10) 

s 

(3-11) 

fl= — S 

such that the functions ipi{n,m,T; P) and ip2{n,m,T; P) satisfy the following system of differential 
equations: 

= (3.12) 

Proof. The proof is standard. For a given s we show how to construct the operator A] i, the other 
cases being similar. 

The derivative of the function ipi{n, m, r; P) with respect to has the following expansions at the 
marked points P^, where by k we denote the appropriate local coordinate kf: 

^jiJi{n,m,T-P) = exp (±X^^>" ) ■ 0(A;^"+") as P ^ P^^, (3.13) 



d 



jMn,m,r;P) = exp ( ^1^-^" ) 0{k'"''-^) as P ^ P^, (3.14) 



,cr=l 



^jiji{n,m,T;P) = exp ( -^r^r ) • 0(A;~™) as P ^ P^ ■ (3.15) 



* \ (T=l 
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Therefore, for an appropriate choice of functions fg^^{n,m,T), the function 

d " 

ij{n, m, t; P) = m, t; P) - ^ fli,^{n, m, r)V'i(n + H, m, r; P) (3.16) 

* fl= — S 

has the following expansions at P^: 

P+, (3.17) 



pr, (3.18) 

and the same expansions fl3.14p - 03.151) at P.^ as ^'ipi{n, m, r; P). Therefore, by (13. 5p this function is 
identically zero on X. Hence, the function ipi{n, m, r; P) satisfies the system of equations (I3.12p . 

Proposition 8 The left ideal of matrix difference operators in DJl that annihilate ip is the principal 
left ideal generated by the operator D. 

Proof. Suppose that A and B are two operators in O that satisfy the following equation: 

A^, + B^2 = 0. (3.19) 

We need to show that there exist elements C,D E D such that A = C{T2 — a) — D(3 and B = 
-CP + D{Ti-a). 

First, we multiply equation (I3.19P on the left by sufficiently high powers of Ti and T2 so that the 
operators A and B become polynomial in Ti and T2. Next, we show that we can eliminate all terms 
containing mixed powers of Ti and T2. Indeed, suppose 

n— 1 

A = aiTlT2~^ + (terms with no T1T2) + (terms of order < n), 
1=1 

n-l 

B = hiTlT2~^ + (terms with no T1T2) + (terms of order < n), 

i=l 

then we can write 

n-l 

A = Y^ [aiTiT^-'^\T2 -a)- hiTiT^-''^0\ + (terms with no T1T2) + (terms of order < n), 

i=l 

n-l 

B = Y^ [biTiT^-'-\Ti -a)- aiTlT^-'-^a] + (terms with no T1T2) + (terms of order < n), 

i=l 

and proceeding in this way, we can eliminate all terms which are not powers of only Ti or T2. Therefore, 
we can assume that A = Ai(Ti) + y42(T2), B = Bi{Ti) + 52(^2), where the Ai, Bi are polynomials in 
only Ti. 



tp{n, m, r; P) = exp ( ± ^ ) • 0(r"^) as P 



^(n, m, r; P) = exp ( ± ^ ) ■ 0(A;-") as P 
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Suppose that Ai = ^"^QOiT^ and Bi = Yl^=o^j'^i- Comparing the singularities in fl3.19p at the 
point Pi, we see that m = n + 1. Subtracting 6„+iT" [(Ti — a)ip2 — Pi'i] from (13.191) . we reduce the 
degree of Bi, and hence of Ai. In this way we can ehminate Ai, and similarly B2. Therefore, we are 
left with showing that ii A = A2{T2) and B = Bi{Ti) are linear polynomials satisfying (13.191) . then 
they can be expressed a.s A = /(T2 — a) — g[3 and B = — f P + g{Ti — a) for some functions / and g, 
which can be easily shown. 

Proposition 9 There exist matrix difference operators BI in 971 such that the following equations are 
satisfied: 

-^D = DA\ + BID (3.20) 

Proof. Equations ([3S]) and fl312D imply that 

d 



dtl 



ip = 0. (3.21) 



Since the operator in the left hand side does not contain derivation in time, it is inside DJt, hence by 
the above proposition it is a left multiple of D, which proves the statement. 

Proposition 10 The equations 

d 

—D + DAI = mod D (3.22) 

dtl " 

define a commuting hierarchy of differential- difference equations. 

We call this system the discrete modified Novikov-Veselov (dmNV) hierarchy. In the next section, 
we give the explicit form of the first two pairs of equations of the dmNV hierarchy. 



4 First and second equations: explicit forms 

In this section, we write down the explicit form of the dmNV hierarchy corresponding to times rl, rf, 
T2 and r|. We give the explicit calculations for tI, the derivations for the other times being similar. 

It is difficult to write down the dmNV as they are defined in (I3.22p . since this involves performing 
division with remainder in a matrix algebra over a non-commutative operator ring. To circumvent this 
difficulty, we notice that the discrete Dirac equation (13.81) . which is a difference equation of degree one 
on the two functions ipi and ip2, is equivalent to a degree two difference equation on one of the ipi or 

^2- 



Proposition 11 Suppose the functions ipi and ip2 satisfy the discrete Dirac equation (EZ^j- Then the 
functions ipi and ip2 satisfy the following discrete Schrddinger equations 



H21P2 



I1I2 - {tia)Ii — 12 + 

T1T2 - {t2a)T2 ^ — Ti + — 



^1 = 

^2=0. 



(4.1) 
(4.2) 
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Proof. This follows from excluding ipi or ip2 from the system (13.81) . 

Conversely, we have an analogue of Proposition 3.4 for the operators Hf. 

Proposition 12 The left ideal of difference operators in O that annihilate ipi is the principal left ideal 
generated by the operator Hi . 

Proof. Suppose that A G D is an operator such that Aipi = 0. Then Proposition 3.4 implies that 
there exist operators C,D E O such that 

A = C{T2-a)-D(3, -Cf3 + D{Ti - a) = 0. 

Expressing C = D{Ti — a){l3)~^ from the second equation and plugging it in to the first, we get that 
A = D{ti(3)^^Hi. The case of il)2 is similar. 

These two propositions allow us to write our hierarchy as a system of rank one difference equations 
of degree two. 

Proposition 13 The discrete modified Novikov- Veselov hierarchy 1^3. 22\} is equivalent to either of the 
following two systems of equations 

-^H, + HiA\ ^ = mod Hi, (4.3) 



-^H2 + H^Ai ^ = Qmod H2. (4.4) 
or! 



We now use this approach to construct the equations corresponding to times r/, rf, and r|. 

The functions and il)2 have the following power series expansions at the marked points 
where by k we denote the appropriate local coordinate kf: 

^i(n,m,r;P) = A;±"exp|^±fjr>'^j ■ (^^it^^{n,m,T)k-'^ e^s P ^ P^, 

^i(n,m,r;P) = /c±-exp ij^r^r • ij2^t(^,m,T)k-''] as P ^ P^^, 



(T=l / \a=0 

00 \ / 00 



(4.5) 



^2(n,m,r;P) = A;±"exp ±^r>M ■ ^t"^""' ""^^^^ as P ^ P^, 

\ a=l / \a=0 J 

^2(n,m,r;P) = A;±-exp|^±fjr2rj ■ (^^xtjn,m,T)k~''^ as P ^ P^, 

where the ^fs{n, m, r) and xfsi''^^ ^1 analytic functions in the variables r, and = 0, xto = 0. 

To make our notation consistent with (I2.8l) - (l2.9p . we denote 

«? = efo, bf = xto (4.6) 

cf = efi, dt = xt (4.7) 

Plugging these expressions into (13. 8p . we see that these coefficients satisfy the following system of 
equations: 
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Ut.±i = <t. + I3xta (4-9) 

tlXta±l=Ktc. + (^xtc. (4.10) 
tlXto.=Kto. + ^xta (4.11) 

Also, since the functions ipi and ?/'2 satisfy the normahzation conditions (]2.18p -( j2.20p . we also have 

a+ar = l, 6J&2" = 1- (4.12) 

We now derive the dmNV equation corresponding to time tI using its equivalent form (14.31) . Let 
/ denote differentiation by tI. We denote A\ ^ = ATi + BT{^ + C and Hi = T1T2 + xTi + 1/T2 + z. 
The equation in time has the form 

- xTi -yT2-z = {T1T2 + xTi + yT2 + z){ATi + fiTi-i + C) mod i/i. (4.13) 

First, we express all of the coefficients of the above equation in terms of the variables af, 6^, a and 
(3. The coefficients x, y, z of Hi were found above in Proposition 4.1: 

. = y=-^, z='f. (4,14) 

To calculate the coefficients of the operator A\^-^^, we use the method of Proposition 3.3. Comparing 
singularities, we see that if 

iiiii t-]^ (1^ (ii 

then the functions ipi ^ind ipi " ^^iV^i ~ BT^^ip are proportional. Hence we can determine the 
third coefficient C = fl^^hy comparing these two functions at either P2 or P2 , which gives us two 
alternative expressions: 

^ .1 1 / dct at , ar 1 I \ 1 / da'^ at ^7 1 \ . n 

We first these expressions by removing the coefficients and c^. From the system (14.8114. 11| we get 
that C2 = (t2^P)(t2^b2) si^nd = — /9/(a6^). Using = the first expression becomes 

^ _ ^2 1/5 , t2 ^6+ at {t2'a){tit2'P) ^ t^^ai t^V/? 



^2-1/5 t2'&2^ t2'f3 at {t2'P){tiH2'a) 

and the second expression becomes 

_ 1 _ /3 a bt af ti(3 t^^af a{t^^(3) 
^ ^^'1'° ^^'a^bf^ tia+ f3{tia) ^ ~af ^ ' 

Expanding the right hand side of (14.131) . we get 

HiAl^ = {tit2A)TlT2 + x{tiA)Tl + [tit2C + y{t2A)] T1T2 + [x{tiC) + zA] Ti+ 
+ [tit2B + y{t2C)] T2 + x(ti5) + zC + y{t2B)T{^T2 + ;z5rf ^ 
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This expression is a Laurent polynomial in Ti and T2 whose terms have degrees i and j in Ti and T2, 
respectively, where i = —1,0, 1,2 and j = 0, 1. We need to express it as a left multiple of Hi plus 
an operator containing terms of degrees (0,0), (0,1) and (1,0). First, to cancel the term containing 
T1T2, we subtract the following left multiple of Hi. 

{tit2A)TiHi = {tit2A)T^T2 + {tit2A){tix)Tl + {tit2A){tiy)TiT2 + {tit2A){tiz)Ti. 

Using (14. ip . (14.151) and (14.81) . we see that the coefficient in front of Tf in this difference vanishes: 

x{tiA) - {tix){tit2A) = -(ti«)^ + {tla)^j^ = 0. 

t^dl t^t2Cii 

Similarly, to cancel the term containing Ti^T2, we subtract 



u D\rr-i -irj yihB) y{t2B) , I y{t2B) , 

vihBjTi y Hi= . T2 + [t^ x) +y(t2B)T^ T2 + [t^ z)T^ , 

h y h y h y 

and using (14.11) . (I4.15p . (14. 8 p and the relation (14.120 . we show that the coefficient in front of 
vanishes: 

Hence, we see that 



^(trM = o. 



HiA\ I = [tit2C + y{t2A) - {tit2A){tiy)] T1T2 + [x{tiC) + zA - {tit2A){tiz)] Ti+ 



+ 



'tit2B + y{t2C)-y^''''^ 



T2 + x{tiB) + zC - ^ {ti\) mod Hi. 

ti y 



ti'y 

Finally, to obtain the evolution equation, we subtract [tit2C + yit2A) — {tit2A){tiy)] Hi from the right 
hand side of the equation, and obtain the following equations: 

-x = x{tiC) +zA- {tit2A){tiz) - X [tit2C + y{t2A) - {tit2A){tiy)] , (4.17) 
-y = tit2B + y{t2C) - 4^ - y [tit2C + y{t2A) - {tit2A){tiy)] , (4.18) 

h y 

- z = xitiB) + zC- y^0p-{ti'x) - [tit2C + y{t2A) - {tit2A){tiy)] . (4.19) 

^1 y 

Since the coefficients x, y, z of H are expressed in terms of a and (3, which are in turn related by 
the equation — = 1, it is sufficient to find one of the derivatives, for example x. Expanding 
the expression for x and using the expressions for the coefficients x, y, z and A, B, C obtained above 
(using the first expression for C in tit2C and using the second one in tiC), we obtain the following 
equation 

% = (4,20) 

which is the first equation of the dmNV hierarchy. 

It seems natural to replace the variables af and 6^ with their logarithms, i.e. to introduce new 
variables af = e"^ and 6^ = e"^. Since a = t2af/a^ = these variables are related by the 

equation 

t2ip - if = hip - ip. (4.21) 
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Writing the evolution equation (14.201) in terms of these new variables, we get 



(4.22) 



To derive the evolution equation for time , we use its equivalent form (14. 4p . The calculations 
in this case are identical to those performed above. In fact, since our problem is symmetric with 
respect to exchanging the marked points and P^, we can obtained the desired equation simply by 
exchanging the functions and 6^ and simultaneously exchanging the shift operators ti and t2 in the 
evolution equation in time (I4.22p . This gives us the following equation: 

In terms of the logarithmic variables, this equation reads 



(4.24) 



The derivation of the equations for times T2 and r| involves similar calculations. For time T2, we 
use the equivalent form (14. 3p : 

- —4 = H^Al^ mod Hi. (4.25) 
(JT2 

Here ^2^1 is a Laurent polynomial in Ti with terms of degree —2 to 3. As above, we successively 
subtract appropriate left multiples of Hi to cancel the terms containing TIT2 for i = 3, —2, 2, —1. At 
every step, the corresponding T{ term vanishes. Finally, canceling the T1T2 term gives us the following 
equation: 

hbt ^ PjtiP) 1 _ PihP) at ^2„+ a| a(tr^/?)(ti/?) t^^a^ 

tia tittj^^^ tia (tiaj^)(tfa]^) (tia)(tfa) tfa^^ tia tiaf ' 

where the functions , b^, , a and (3 in the equation satisfy the following relations: 

a = ^ = ^, a'-(3' = l, hct = act + f3{ti'mi'cit). (4.27) 



5 Explicit formulas 



In this section we give explicit formulas for the functions tpi{n, m, r; P) in terms of the theta-functions 
of the surface X. Choose a basis aj, bj, j = 1, . . . , g of Hi{X, Z) with canonical intersection form, i.e. 
such that ttj o = 0, bj obk = 0, a, 06^ = Sjk- Let B be the period matrix of the curve X with respect 

r1 ifl-OT-Cin-f- 1 n 1 o ^X^tn ^\^Tvr\ l^ini^ "ciri-i-Vi T-\/~\lciCi o m J ^ 



to this basis. Let and denote Abelian differentials of the third kind with poles at Pj^ and P^- 



Q] = d{kfr' {Tkf + 0(1)) as P ^ 

which are normalized to have zero periods over the a-cycles. Let VL^ denote Abelian differentials of the 
second kind with poles at P^ and principal parts 

= d{kf)-' {Ts{ktr' + 0(1)) as P ^ Pt, 
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and with zero a-periods, and which are odd with respect to the involution a. It is a standard fact that 
these differentials exist and are unique. Let U} and denote the vectors of the 6-periods of these 
differentials: 



SI'. 



Choose a base point Pq E X away from the marked points and the divisor D, and let A: X ^ J{X) 
denote the Abel map with base point Pq, where J{X) is the Jacobian variety of X. Let d{z\B) denote 
the theta function of J{X) for z G C^. Introduce the functions 



n{P) 



r2{P) 



ejAjp) - A(p+) - MPd - K\B)e{A{p) - YJitl a{p.) + a(p+) - k\b) 

0{A{P) - ELi m) - K\B)9{A{P) - Eg A{P,) - K\B) 
e{A{P) - A{P^) - ZI2 A{P.) - K\B)e{A{P) - Eg A{P.:) + A{P^ 



K\B) 



d{A{P) - Eg A{Pd - K\B)d{A{P) - Eg A{P,) -K\B) 

By construction, these are meromorphic functions on X whose pole divisor is D = Eg^ Pi ^^^1 whose 
zero divisors are P2 + -Di and P^ + D2-, respectively, where Di and D2 are some divisors of degree g. 
We define the functions ipi and il)2 by the following formulas: 

00 2 



n 



m, r; P) = ri{P)Ci{n, m, T)Fi{n, m, r; P) exp 

where the function F{n, m, r; P) is defined as 

/ 00 2 

e A{P) - + nUl + mUl + 



Po 



(5.1) 



Fi{n,m,T;P) 



s=l i=l 



e{A{P)-A{D,)-K) 



and the path of integration in the exponent is the same as in the Abel map in Pj. By construction, 
these are single-valued functions on the curve X, having the required meromorphic and exponential 
singularities at the marked points, and having pole divisor D away from the marked points. 

The constants Ci{n,m,T) are determined by the normalization conditions fl2.18p - fl2.20l) . Choose 
paths of integration 7^ : [0, 1] ^ X from Pq to P^ and a path 7 from Pq to cr(Po). We assume that the 
integration path in ipi{P) is 7^ and that the path in ipi^a^P)) is 7 followed by the image of 7^ under a. 
Writing out the expression for ipi{P)ipi{a{P)) using (15. ip . we see that we need to choose the constants 
Cj(n, m, r) as follows: 



Ci{n,m,Ty 



ri{Pi^)ri{Pi )Fi{n,m,T]P^)Fi{n,m,T]Pi ) exp 



00 2 „ ' 

s=l i=l J-1 



(5.2) 



where the path of integration in the Pi(n, m, r; P^ ) factor is 7 followed by o"(7j), and the constants l\ 
and Pf are the principal values of the integrals of l^j^ and fig along the path — 7j + 7 + cr(7j): 




<7{7(t)) 



1,2. 



(5.3) 



<7(7(0)) 



Finally, we choose the signs of Cj(n, m, r) in such a way that the functions z/'i satisfy the equation 
fIl2QD. 
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